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We present theory and experiments for the force-distance curve F{zo) of an atomic force mi- 
croscope (AFM) tip (radius i?) indenting a supported fluid bilayer (thickness 2d). For reahstic 
conditions the force is dominated by the area compressibility modulus ka of the bilayer, and, to an 
excellent approximation, given by = 'jTK,ARzo/{2d — zq)'^. The experimental AFM force curves 
from coexisting liquid ordered and liquid disordered domains in 3-component lipid bilayers are well- 
described by our model, and provides ka in agreement with literature values. The liquid ordered 
phase has a yield-like response that we model as due to the breaking of hydrogen bonds. 

PACS numbers: 87.16. D-, 87.16.dm, 87.80.Ek, 68.37.Ps 



I. INTRODUCTION 

Atomic Force Microscopy (AFM) [l| has become a 
standard tool for imaging surfaces at high resolution and 
probing local mechanical properties [Ij. Force-distance 
curves for indentation of AFM tips have been used to 
characterize the mechanical properties of biological mem- 
branes and the usual approach is to approximate 
the bilayer as an elastic solid undergoing a Hertzian con- 
tact @-@|- However, at physiological conditions most bio- 
logical membranes are in a fluid bilayer phase [H , whose 
free energy is described by a bending modulus k and 
the area compressibility modulus ka- These are exper- 
imentally accessible th roug h, for example, micropipette 
aspiration experiments [ifll , which give the average value 
of the elastic moduli over the whole vesicle. However, 
biological membranes often have different local composi- 
tions, and thus different local mechanical properties and 
physiological functions. 

Despite the growing use of AFM to study lipid bilay- 
ers, the flexibility of using it to measure local mechanical 
properties has not been fully exploited. An AFM tip 
can bend a freely suspended membrane, and compress 
a supported membrane. In recent work, Steltenkamp et 
al. [ll| showed how to extract the bending modulus of 
lipid bilayers from AFM force-distance curves for bilay- 
ers deposited over well defined sized holes (indentation of 
'nanodrums'), in which they could safely ignore area com- 
pression due to the lack of a supported surface. Another 
issue neglected in previous AFM studies is the double 
leaflet form of lipid bilayers, which is known to influence 
the dynamics of fluctuations [l^]- Since an AFM tip in- 
duces an asymmetric response in a supported bilayer [isj , 
the distinction between the two leaflets will be important 
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to accurately model the mechanical response. 

In this paper we consider the force-distance curves ob- 
tained by indenting an AFM tip into a fluid bilayer sup- 
ported on a solid substrate. The force-distance curves 
arc calculated from a static analysis of the deformation 
of the two leaflets and differs from usual Hertzian result 
of the deformation of elastic bodies. We analyze exper- 
iments on a dioleoyl-phosphatidylcholine (DOPC) - egg 
sphingomyelin (SM) - cholesterol (CHOL) phase sepa- 
rated supported bilayer, which is a model mixture repre- 
sentative of typical in vivo membranes [l^ |. For certain 
composition ratios of the components, this system spon- 
taneously phase separates into coexisting liquid ordered 
(Lo, rich in SM and relatively thick because of strong ne- 
matic order in the acyl tails) and liquid disordered (L^, 
rich is DOPC and relatively thin because of the more 
disordered tails) phases. We show how to determine the 
area compressibility moduli of the coexisting Lq and 
phases of a single sample, and find values in agreement 
with literature values. To the best of our knowledge, this 
is the first time that the area compressibilities of the two 
coexisting compositions in fluid bilayers have been ex- 
tracted directly. This technique should prove invaluable 
for studying the composition dependence of mechanical 
properties in lipid bilayers, and can be easily extended 
to consider more complex interactions between AFM tip 
and the bilayer. 



II. THEORY 

We consider a supported fluid lipid bilayer of thickness 
2d probed by an AFM tip of radius R in contact mode, 
which measures the force as a function of the depth zq 
from the unperturbed surface of the layer (see Fig.[lja)). 
At typical AFM speeds the viscous forces are negligible. 
We begin by assuming a hard contact interaction between 
the tip and the membrane. Electrostatic repulsion from 
the charged double layers and van der Waals attraction 
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are included later in the paper when comparisons are 
made with experiments. Since the lipid bilayer is not 
anchored it remains tension free. We assume that the 
volume is conserved at the molecular level: as the tip 
penetrates the bilayer it occupies a volume SV, so that 
SN = 5V/ (aod) lipids are expelled into the surrounding 
bilayer. Here, og is the area per lipid in the absence of the 
AFM tip. The surface area increases by SA, due to the 
curved spherical surface of the AFM tip, and the increase 
in area per head group da = a — is given by 



A + SA JA + SA 

da = — — ao ~ aod 



N-SN 



V-SV 



(1) 



This increased area induces an elastic cost due to the 
stretching elasticity of the lipid leaflets. We calculate 
this not by averaging over the entire spherical cap, but by 
considering small increases in radius dr, and evaluating 
Sa/ao at each r (Fig. [1]). 





FIG. 1: (color online) (a) Schematic geometry of an AFM 
tip of radius R indenting a fluid bilayer of thickness 2d by an 
amount zq. The leaflet dividing surface at hb{r) is shown as 
a solid line, (b) Cylindrical volume elements of depth z{r) at 
the distance r from the center of the tip. (c) The area 27rr ds 
of such an element in contact with the lipid. 

We assume that both leaflets have the same area per 
lipid oo and stretching modulus ka/2 for lipid head 
groups on a flat surface. This should be valid in the 
absence of specific interaction of the lipid with the sub- 
strate, although experiments have shown that the sur- 
face often does have specific interactions [l^. The head 
groups in the top leaflet are forced to lie on a curved sur- 
face below the AFM tip. This affects both the area/lipid, 
and the stretching modulus for the lipids on the top 
leaflet. We model the local lipid free energy as a sum of 
a surface energy and a harmonic tail stretching, ^(i) ~ 
70^/ cosd + a/aj^, where 7 is a surface tension, a penal- 
izes tail stretching, and a^ is the projected area for leaflet 



thickness L. Here, 9{r) = sm~^ (r/R) (Fig.[T]:) is the tilt 
angle of the lipid surface. Minimizing g at fixed lipid 
volume 1; ~ LaL leads to an effective stretching modu- 
lus ~ sec^^^ 9 and an effective area per lipid 
flg = ao cos^/"^ for the top leaflet. 

Using the modiflcations due to the curved surface for 
the head-groups in the top leaflet, the excess free energy 
due to the increase in the area per lipid during indenta- 
tion is 



d^r 



sec2/3 9{r) I — 



6a 



Sa 

ao/b 



(2) 



where t and b refers to the top and bottom leaflets and the 
integration extends over both leaflets. The lower leaflet 
will generally deform to accommodate the large energy 
change due to removing too many lipids from the upper 
leaflet. We let the lower and upper leaflets have thick- 
nesses hh{r) and ht{r) respectively, with hb{r) + ht{r) — 
2d— zo + R[l — cos 6{r)] . The area changes at each radius 
r are given by (Fig. [T]) 



/ t 
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a, 
da 



dsec^/^9{r) 



htir) 



hb(r) 



(3a) 
(3b) 



The measure is d^r = rd(j)dr. The dividing surface hb(r) 
is determined by minimizing the free energy at each r. 
For equal stretching moduli in both leaflets, an explicit 
solution for hb{r) is possible for small tilt angle. 
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2d- zq + R {1- cos 9{r)) 
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(4) 



For realistic values zq ~ 2nm, d ~ 3nm and R ^ lOnm, 
this approximation introduces less than 0.1% error in 
hb{r) for the entire range of r. We use this approximation 
in the rest of the paper to derive analytic expressions for 
the free energy and force. 

Using hm{r) from Eq. 21 the free energy is 



2G(zo) 
ttkaR^ 



-2/3 



,1/3) 



-\ 2 



r-1/3 



[2-lf + |(l-^)] 



dx. 



(5) 



and the force on the AFM tip is given by F = dG/dzQ. 
We use the numerical force derived from Eq. [5]when per- 
forming flts to the data. For small penetrations zq the 
force can be written as 



F 



ttkaR 



d 

3R 



(4)l(7r-°<^.')-(«) 



A surprisingly simple function that fits the entire exper- 
imental range of forces, correct to within a few percent 
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for R = 3d and much better for larger R, is 

/ \ 2 



4 \2d- zo 



(7) 



The force diverges as zq approaches 2d because the 
area/Upid diverges in order to preserve molecular vol- 
ume. The quadratic free energy (Eq. [5]) is no longer 
valid there. Experimentally this divergence is preempted 
by pore formation (see below). 

For comparison, the contact force between two solid 
(elastic) bodies much larger than the radius of contact 



(Hertzian contact) is F 



3/2 



1151. More relevant for 



AFM experiments, the force to indent a finite elastic 
layer scales as F R^z^/d^ if the layer is bonded to 
the substrate and F ~ Rz^/d if the layer can slip 0. 
The response of fluid bilayers in eqn. [S] scales diff'erently 
than all of these scenarios, and force is proportional to 
the area compressibility modulus instead of the Young's 
modulus. For realistic experimental values the region of 
validity for this quadratic behavior is limited, as shown 
in Fig.H^b). 
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FIG. 2: Scaled force F / (jikaR) as a function of tip depth 
Z(i/d (solid line) for tip radius R/d = 3, according to Eq. [S] 
The dashed line shows only the leading quadratic term in 
zo/d, Eq. [S] Inset: Behaviour at small zo/d. On this scale 
the approximation of Eq. ((7|) is indistinguishable from the 
numerical solution of Eq. ([5]). 



III. EXPERIMENTS 

To test the theory we performed experiments on a sup- 
ported bilayer comprising DOPC, SM and CHOL at over- 
all molar ratios DOPC:SM:CHOL=40:40:20. At room 
temperature this system phase separates into coexisting 
DOPC-rich liquid disordered and SM-rich liquid ordered 
domains (Fig. The hydrocarbon tails have large ne- 
matic order in the Lq phase, leading to a thicker bilayer 



and higher area compressibility modulus. In contrast, 
the tails have lower nematic order in the L^; phase with 
concomitant smaller thickness and lower moduli. DOPC, 
CHOL (purchased from Sigma) and egg SM (purchased 
from Avanti) were dissolved in chloroform, dried under a 
stream of argon for 30 minutes, and then vacuum desic- 
cated for 30 minutes. The lipid was rcsuspcnded in PBS 
buffer at pH 7.4 to a concentration of 1 mg/ml by vor- 
tcxing. To make small single unilamellar vesicles (SUVs), 
the cloudy lipid suspension was tip sonicated (IKA, U50) 
at less than 5° C for 25 mins (until the solution became 
clear). The mica (Agar Scientific Ltd.) surface was in- 
cubated with the SUVs at at 50° C and cooled down to 
room temperature in a incubator over 15 minutes. Af- 
ter Ih, the sample was gently rinsed with PBS buffer to 
remove any excess vesicles. 

Force measurements were performed at 27° C in PBS 
buffer using a Nanoscope IV Mulitmode AFM (Veeco) 
equipped with a temperature control stage, using can- 
tilevers (NP, Veeco) with nominal spring constants of 0.12 
N/m. Spring constants were measured using the thermal 
noise method [l^ in air, and optical lever sensitivity de- 
termined against a clean mica surface. The force curves 
analysed in this paper were all taken from a single force- 
volume map of the phase separated bilayer shown in Fig- 
ure. 21 and exported using Nanscope software v5.12r30. 

Scanning electron microscopy (Camscan series III, 
FEG-SEM operating at 5 kV with magnification 160k) 
was used to measure the tip radius. Inset of Fig. 2] shows 
the tip image with dashed lines along the edges of the 
four pyramidal faces. The end of the tip can be approx- 
imated as spherical. The drawn circle (Fig. |4] inset) has 
a radius of 10 nm. The contrast of the image is poor. 
Consequently the uncertainty of the exact value of the 
radius is large. In our analysis we consider the tip radius 
to be i? = 10 ± 5 nm. 

Fig.[3]shows a tapping mode image of the bilayer along 
with a one dimensional cross section. There is a ^ 5 nm 
thick Lrf matrix enclosing ^ 6 nm thick Lq domains (The 
heights reported here include the thickness of any wa- 
ter layer between the bilayer and the mica surface). The 
composition of the two phases were determined by fol- 
lowing the tie lines on the ternary phase diagrarn, which 
were calculated using Atomic Force Microscopy |17| : 



Ld : (DOPC:SM:CHOL) = (68 : 27 : 5) 
Lo : (DOPC:SM:CHOL) = (3 : 71 : 26) 



(8b) 



Phase diagrams on similar ternary mixtures have been 
calculated using NMR, and the compositions of the 
liquid-disordered and liquid ordered phases are similar 
[18| . The uncertainty in the compositions from place- 
ment of the tie lines is estimated to be less than 2% of 
the quoted values. For both Lo and Ld phases, force 
curves from contact mode AFM were used from at least 
10 different measurements from different points within 
different regions ('patches', as in Fig.[3|) of the sample. 

Fig. 0] (symbols) shows the force curve for the DOPC 
rich bilayer in the liquid disordered phase. Besides the 
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FIG. 3: (color online) Phase separated lipid bilayer with liquid 
ordered and liquid disordered domains, (a) Tapping mode 
AFM image showing the height profile of the bilayer. (b) One 
dimensional section along the dashed line in (a). 




FIG. 4: Force-distance curve for a DOPC rich bilayer in the 
liquid disordered phase. The data points are from AFM ex- 
periments and the line is a fit for the theoretical prediction 
with ka = 0.12 N/m. Inset: SEM image of the tip to measure 
the tip radius. The dashed lines are along the pyramidal face 
edges. The circle drawn at the end of the tip has a radius 
10 nm. 



stretching contribution considered so far, the tip experi- 
ences an attractive force due to van der Waals interaction 
and a short range repulsive force due to the electric dou- 
ble layers on the tip surface and the membrane top sur- 
face. In principle the van der Waals interaction can be 
calculated from a knowledge of the dielectric constants 
of the tip, membrane, and the PBS buffer [l9| . Simi- 
larly, the repulsive interaction can be estimated by know- 
ing the detailed charge distribution and solving Poisson- 
Boltzmann equation. Phenomenologically, we model the 
van der Waals attraction as an interaction energy of the 
form —A/S^^ between volume elements of the tip and the 
membrane separated by a distance ^. 

Since these forces are short-ranged, we consider the tip 
as a sphere and for the volume integration over the mem- 
brane, consider the membrane as infinitely thick. We 
further assume that the repulsive interaction is strong 
enough to avoid adsorption. As the tip approaches the 
bilayer, the bilayer deforms. The extent of the deforma- 
tion is governed by the minimum of the stretching free 
energy and the long range interactions (van der Waals 
and screened Coulomb). We assume that the deforma- 
tion can be modeled as hard interaction from a tip with 
radius Rc larger than the physical tip radius R. Although 
A is poorly known and depends on the detailed dielec- 
tric properties of the membrane, its precise value only 
changes Rc and controls the details of the force near 
contact. For deeper contact the force is overwhelmingly 
dominated by the stretching modulus ka, so that the 
force-distance curves yield the same ka, independent of 
A. 

The drawn line in Fig. |4] shows the fit from our the- 
oretical analysis, using a downhill simplex method (20j 
to minimize the mean square fractional deviation of the 
prediction from the experimental data over the fitted 
range. The best fit for the compressibility modulus for 
R = 10 nm is ka ~ 0.12 N/m. Because of the uncer- 
tainty in the tip radius, the range of ka for R between 
5nm and 10 nm is between 0.25 N/m and 0.08 N/m, re- 
spectively. Our estimate compares well with the litera- 
ture values Ka _^0.13 — 0.6 N/m from osmotic pressure 



measurements 



21| and ka = 0.18 ± 0.04 N/m from mi- 
cropipette aspiration of GUVs [l^] made of pure DOPC. 
Our model provides an excellent fit until 2(1— zq ~ 2.5nm, 
at which point the elastic energy of the deformed bilayer 
overcomes the cost of forming a hole [1^ and the tip 
abruptly penetrates the full bilayer. 



IV. RESPONSE OF LIQUID ORDERED 
DOMAINS 



The AFM force curves for the SM rich liquid ordered 
phase are qualitatively different from those in the coex- 
isting liquid disorded phase (Fig. [5]). The initial deforma- 
tion (5 nm < 2d — zq < 6 nm) shows a high modulus con- 
sistent with the tightly packed character of the Lq phase. 
Around 2(i — zq ~ 5 nm the response shows a crossover to 
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FIG. 5: (color online) Force-distance curve from AFM (sym- 
bols) in the SM rich liquid ordered phase superposed with 
two separate theoretical fits (lines, using Eq. [SJ involving two 
different ka at small and large tip penetrations zq. Inset; 
Closeup of the crossover region. 



a much lower modulus. The symbols in Fig. [5] are from 
12 separate force-distance measurements. While the ex- 
perimental data fall on the same curves away from the 
crossover region, the transition from stiff to soft behav- 
ior occurs at different values of zq, which may be due to 
either the stochastic behaviour of an activated event or 
flucuations in composition from region to region. 

We first attempt to model these force curves as due to 
an effective stretching modulus that differs for small and 
large penetrations far from the crossover region. Hence 
we fit the data at small penetration {2d — zq > 5.3 nm) 
and large penetration {2d — zq < 3.5 nm), with effective 
stretching moduh according to Eq. ([5]). For R = 10 nm, 
the small zq fit gives ka = 1.1 N/m. Recent experi- 
ments on a bovine brain SM and CHOL equimolar mix- 
ture found KA = 2.1 ± 0.2 N/m [H. Since egg SM (16:0 
SM) has shorter fatty acid chains than does bovine SM 
(18:0 SM), and the current composition has compara- 
tively smaller amounts of CHOL, we expect the mem- 
brane to be softer (smaller ka), as found. 

The large penetration {zq) region has a stretching mod- 
ulus KA = 0.05 N/m, which is much closer to that of the 
Ld phase shown in Fig.|3]than the unperturbed Lq phase. 
The AFM tip forces the bilayer immediately below it to 
decrease in thickness, which thus destroys the strong ne- 
matic order of the Lq phase and induces a yielding or 
phase transition of the Lq phase into an phase. It is 
likely that the composition of this induced Ijd phase dif- 
fers from that of the phase that characterizes equilib- 
rium coexistence far from the AFM tip (Eq. I5al) , because 
of slow kinetics of composition changes under the AFM 




2 d - Zq (nm) 

FIG. 6: (color online) Force-distance curve from AFM (sym- 
bols) in the SM rich liquid ordered phase superposed with a 
microscopically-motivated fit that accounts for a separate en- 
ergetic contribution from hydrogen-bond breaking (solid line, 
based on Eq. [9]). Also shown are the separate contributions 
from the van der Waals interaction (dotted line) , from the hy- 
drogen bonds (dot-dashed line) and the area compressibility 
term (dashed line). 



tip. Our separate fits to extract ka suffer from narrow 
available fit window (~ 0.3 nm) for small zq and the lack 
of small force data for the large zq fit. Also, this proce- 
dure does not elucidate the reason for two distinct elastic 
regions separated by a crossover. 

To understand the qualitatively difference force re- 
sponses of the Lo and phases, we propose a 
microscopically-motivated model. SM has both hydro- 
gen bond donor and acceptor gro ups and is known to 
form inter-SM hydrogen bonds |25l . |26| . The free energy 
in the phase, as represented in Eq. [21 is dominated 
by solvent and tail packing entropies. Hence, to describe 
the Lo phase we separately include the short range en- 
ergy of hydrogen bond breaking through a sim|)le Morse 
potential: U{b) ~ Ed [1 — exp {—{b — &o)/Am)] , where b 
is the separation between the donor and acceptor group 
and 6o is the equilibrium separation. For typical hydro- 
gen bonds the dissociation energy Eu ^2 — 7KCal/mol 
and the range A„ 0.02 - 0.07 nm [13, [111. For smaU 
changes in area/lipid and affine deformation the contri- 
bution to the free energy from distortion of hydrogen 
bonds is approximately 



h^{r) 



(9) 



Here, Chb is inter-lipid hydrogen bond dissociation energy 
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per area and A ~ 2>^^^\ra/ y/oo for hexagonal arrangement 
of the Upids. As before, the total free energy, now com- 
prising contributions from Eqs. (0 IH]), is minimized at 
each r to find the dividing surface between the leaflets 
and the force is calculated from F ~ dG/dzQ. 

In the limit of small penetrations this model gives an 
effective stretching modulus k'^ in the Lq phase of 

nf^n^ + i^, (10) 

where ka is thus the stretching modulus of the hd phase 
that is left after the Lq phase has been destabilized and 
there is no remaining hydrogen bond contribution. The 
fit to the data is shown in Fig. [51 The stochastic nature of 
the force curves near the rupture point {2d — zq ~ 5.7), 
limits the ability to obtain excellent fits. Our fit gives 
KA = 0.13 N/m, Chb = 0.006 N/m and A = 0.1. 

Assuming an area per lipid gq ~ 0.6 nm^, the fitted 
value A implies that the range of the Morse potential is 
Am ~ 0.06 nm. Simulations show about 0.5 hydrogen 
bonds per lipid in SM bilayer [25| . Assuming an average 
hydrogen bond energy of 3.5KCal/mol, our value for ens 
give 0.4 hydrogen bonds broken per lipid. The initial 
deformation is dominated by the contribution from the 
hydrogen bonds, and the corresponding force curve leads 
to an area compressibility modulus k'^J'^ ~ 2.7 N/m. 

V. DISCUSSION 

We have assumed a static force response, despite typ- 
ical tip velocities vup — lO^nm/s. We can estimate the 
correction due to finite tip velocity by considering the 
dissipation from two dimensional viscosity r] of the lipid 
layer. The dissipative force is found to be 

i]TTZo{2R- zo) 
Fd{zo)^ 2d? 

The two dimensional shear viscosity for fluid bilayers is 
expected to be of the order of 10^^°N-s/m [2^, leading to 
F]j ~ 10~^ nN, much smaller than the elastic contribu- 
tions. Hence our static approach is sufficient to describe 
the AFM force-distance curves on fluid lipid layers. 

In our calculations we have assumed that the two 
leaflets have the same area compressibility and preferred 
area per head group. This may not always be the case, 
because of surface interactions [isj : for example, sup- 
ported bilayers often have different melting temperatures 
than their counterparts in giant unilamellar vesicles. In- 
corporation of asymmetric membranes into the model is 
straightforward, although more complex. We have also 



neglected splay or bending energies. Part of the elastic 
cost of this is already included in the increased area/lipid 
against the curved surface, in Eq. [21 but there may also 
be an additional negligible free energy cost due to the 
splay of the lipid tails, through the bending modulus of 
each leaflet. 

In our analysis, the initial deformation for the Lq phase 
is described in terms of the stretching of hydrogen bonds. 
This can be explicitly tested by performing experiments 
with varying concentrations of SM or using chemicals 
that disrupt hydrogen bonds. However, this is beyond 
the scope of the present work. 

Evidently, local applied pressure can melt the liquid 
ordered phase into the thinner phase, which is not 
surprising. We have proposed an explicit microscopic 
mechanism in terms of breaking hydrogen bonds that are 
implicated in stabilizing Lq phase. An alternative and 
more general description could include a Landau theory 
for the free energy of the Lq-L^ phase transition, with 
local pressure p added as an external field to destabilize 
the Lo phase, AG ~ pij:, where V' is an order parameter 
proportional to thickness whose value decreases upon a 
transformation to the hd phase [s^l- The phase trans- 
formation would then occur first at constant composi- 
tion, and then one may expect the composition to change 
slowly as the external force changes the local preference 
for the different lipid species. The subject of kinetics and 
composition as a function of applied pressure is interest- 
ing and important, and we leave this for further work. 

In summary, we have presented, and validated by ex- 
periments, a theory for describing the force distance 
F(zo) relationship for AFM experiments on fluid bilayers, 
which leads to a remarkably simple expression for i^(zo), 
Eq. ([7])- This provides a method for flnding the area 
compressibility modulus and the amount of inter-lipid 
hydrogen bonds of fluid bilayers. The agreement with 
the existing literature values for the area compressibility 
is excellent. The main uncertainty in our prediction is 
due to the uncertainty in the tip radius R. However, the 
simple linear dependence on R means that relative mea- 
surements taken with the same tip can be compared very 
accurately. 
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